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We present the analysis on the quadrupole phases of the Internal Linear Combination 
map, ILC(I) and (III) derived by the WMAP team (1 and 3-year data release). This 
approach allows us to see the global trend of non-Gaussianity of the quadrupoles for 
the ILC(III) map through phase correlations with the foregrounds. Significant phase 
correlations is found in between the ILC(III) quadrupole and the WMAP foregrounds 
phases for K-W band: the phases of the ILC(III) quadrupole £2,11 £2,2 and those of 
the foregrounds at K-W bands $2,1, $2,2 display significant symmetry : §2,1 + < I > 2,1 — 
€2,2+^2,21 which is a strong indication that the morphology of the ILC(III) quadrupole is 
mere reflection of that the foreground quadrupole through coupling. To clarify this issue 
we exploit the symmetry of the CMB power, which is invariant under permutation of the 
index m = 1 <^ 2. By simple rotation of the ILC(III) phases with the same angle we reach 
the phases of foreground quadrupole. We discuss possible sources of phase correlation 
and come to the conclusion that the phases of the ILC(III) quadrupole reflect most likely 
systematic effects such as changing of the gain factor for the 3-year data release with 
respect to the 1-year, rather than manifestation of the primordial non-Gaussianity. 
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1. Introduction 

After the release of the Wilkinson Microwave Anisotropy Probe (WMAP ) 1-year 
results (hereafter WMAP I) Refs. 1, 2, 3, 4 the issue of non-Gaussianity of the CMB 
has attracted great attention. In the papers Refs. 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 
various kinds of methods have been employed and departure of Gaussianity have 
been detected. These features can be of primordial origin (Refs. 8, 16) or they could 
be related to foreground residuals (Refs. 17, 18, 19) or any sort of systematic effects 
(9). 

Recently the WMAP team released the 3-year data (hereafter WMAP III) 
Ref. 20 producing Internal Linear Combination map (LLC (III)) available for sci- 
entific analysis at the multipole range I < 10. The WMAP team performed the 
analysis of Gaussianity for the ILC(III) signal that does not include harmonics 
£ < 10 and showed that the CMB signal is Gaussian Ref. 21. However, as it was 
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cautioned by the WMAP team, the alignment and planarity at multipoles £ = 2,3 
and 5 mentioned in Refs. 8, 9, 11, 12, 22 still present in the map (see recent inves- 
tigation in Rcf. 23). 

If breaking of statistical isotropy in the low multipoles is of primordial origin, a 
very fundamental issue is consequently raised: what cosmological model can provide 
such a peculiar structure in the CMB angular distribution on the sky ? In Refs. 8, 
24 and 25 it has been pointed out that breaking of statistical isotropy in the low 
multipoles can be explained and in favor of the Bianchi Vllh model. However, 
Helling, Shupp and Tesilleanu (Ref. 26) have also pointed out that I — 2, 3 and 
5 have anti-alignment in the direction of non-cosmological dipole and argued that 
accurate technique of the CMB signal subtraction probably needs to be non-linear 
(with respect to dipole treatment) to prevent any residuals from the dipole in the 
map. 

In addition, methods of foreground separation for the CMB signal very likely 
leave foreground residuals in the CMB map, a problem that can be illustrated by 
cross-correlation between the derived signal and the foregrounds. For the WMAP I 
data this problem was already investigated in Refs. 27, 28, 17, 18, 29, 30. 

For WMAP III data, it has been shown in Refs. 31, 18 that the low multipoles in 
ILC(III) clearly display significant cross-correlation with the derived foregrounds, 
while Ref. 32 showed that negative peak at b = —57°, I = 209°) has no analogue in 
the foreground map. In Refs. 33, 10 analyzing the WMAP map-making algorithm 
(MMA) it has been pointed out that it could produce some residuals from the non- 
cosmological dipole (related with the motion of our Galaxy in the Local Group) 
and any non-Gaussian features can "naturally" arise as the result of the MMA. 
Summarizing present status of non-Gaussianity and statistical anisotropy of the 
CMB for low multipoles I < 10, we may say that the origin of these peculiarities 
still is uncertain. 

In this paper we go back to the question raised in Refs. 13 and 34, if the 
WMAP quadrupole is cosmological ? To answer this question, we will implement 
the analysis of the LLC (I), LLC (III) and de Oliveira-Costa & Tcgmark (Ref. 35) 
quadrupole and compare their phases with those of the WMAP III foregrounds 21 : 
combined maps from the synchrotron, free-free, dust emission at K-W band. We will 
show that combination of the ILC(III) phases £2,™ and the foreground phases <I>2,m 
follow the intriguing relation £ 2 ,m + 3>2,m — const with accuracy, e.g. within 0.014 
(rad) for the foreground at K band. Since phases of the signal is closely related to 
morphology of the maps (Ref. 36), detected correlation between the LLC (III) and 
foreground phases allow us to conclude that it is most likely due to systematic ef- 
fects (such as estimation of the gain factor and foreground separation) rather then 
primordial one. 



http: //lambda. gsfc .nasa.gov/product/map/m_products . cfm 
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2. Phase analysis of the WMAP data 

For the statistical characterization ol the CMB temperature anisotropies on a sphere 
where 6 and ip are polar and azimuthal angle of the polar system of coordinate, it 
is useful to express AT '(6, p) in terms of spherical harmonics Ye m (8, ip): 

£max I 

AT(6>, p) = y~] ^2 \ a e,m\exp(i(f>i tm )Y( m (e,ip) (1) 
e=o m=-e 

where \ae tT n\ and (f)i^ m are the amplitudes and phases of the coefficients ae tT n, respec- 
tively, and |m| < I. AT (9, <p) can be full-sky signal, such as the WMAP frequency 
maps from K to W band. We denote for the ILC and the foreground coeffi- 
cients as ct.rn = |Q, m |exp(i& >m ) and F^ m = \F^\ exp(i<f>^J, respectively. Here 

and are the corresponding phases and the index j = 1 — 5 marks the 

WMAP frequency band for K, Ka, Q, V and W bands, respectively. In standard 
cosmological models (i.e. those involving the simplest forms of inflation) these tem- 
perature fluctuations constitute a realization of a statistically homogeneous and 
isotropic Gaussian random field. The statistical properties of Gaussian random fields 
are completely determined by the power spectrum 

(«£m «£',„') =C t Su' S mm'' ( 2 ) 

where the angle brackets indicate ensemble average. For any single realization of 
signal, the angular power spectrum is estimated 

m— — £ 

As the signal is always real, the conjugate properties of the spherical harmonic 
coefficients allow us to write down the angular power spectrum as 

^)~Ko| 2 + ^Ek m | 2 - (4) 



2.1. Linear phase correlation method. 

To investigate different correlation of phases for different signals (the CMB and 
foregrounds) we need to implement as much as possible methods to detect corre- 
sponding correlations and then to show their possible sources. One of the possible 
method is to draw different linear combination of the phases in form 

G,, m = 5>v£ (5) 

3 

where ctj is (—1, 0, 1) and ipj ^ are the phases of different signals j. The basis of this 
method has simple motivation from the analysis of the non-Gaussian random process 
recently performed in Refs. 37 and 38. If non-linearity of the random process 5(x), 
which leads to its non-Gaussianity has simple quadratic form, the standard method 
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of the correlation analysis is the bispectrum. It is defined as Fourier transform of 
the tree-point correlation function 

C(n,r 2 ) = (S(x)S(x + f{)5(x + r$)> (6) 

Generally speaking for any non-Gaussian processes we need to know all high order 
cross-correlations or polyspectra, which are negligible for the Gaussian random pro- 
cess. However, we may characterize this sort of Non-Gaussianity in terms of phase 
correlation. For process characterized by the bispectrum only, corresponding equa- 
tion for the phase coupling, which maximizes the bispectrum, has a form (Refs. 37, 
38) 

G(fc) = 2% - e 2V G(k,p) = e~ k + e p e p+% (?) 

where 9^ is the Fourier phase of the random process S(x). 

As one can see from Eq.(7), linear combination of the phases characterizes non- 
linear coupling of the modes tested by bispectrum or in general case by polyspec- 
trum. This is the reason why we introduce the functional Gi. m to investigate dif- 
ferent high order moments of the CMB-foreground cross-correlation. 

2.2. Symmetry of the power spectrum. 

Let us discuss the properties of the second term in Eq.(4). This part of the 
power spectrum is invariant under the following transformation d^ m = pc e m ' , 
de, m =o — Q im= o, where p is the permutation operator for m < I. For the 
quadrupole component 

'd 2 ,A = ( 1 \ fc 2 ,i 
d % i) \ 1 J \c 2 , 2/ 

i.e. this operator simply exchanges m — 1 and 2 without altering the power spectrum 
C{tj. However, the map after such implementation of the operator p 

AT(e,ip) = J2 \d e , m \e**^Y im (0,v) (9) 

£=2 rn=-i 

is not invariant. As one can see from Eq.(8), the phases of the cfe,™, and ci. m coef- 
ficients correspond to the following equations: ^24 = £2,2 and ^2,2 = £24- What is 
important is that if cg^ m constitute a Gaussian random field the phases should 
have no correlation with the foreground phases $£ jm , which should also be the case 
for the ^i t m as well. In Fig.l we plot the ILC(I) quadrupole in galactic coordinates 
and the corresponding transformed map by the permutation operator p, which have 
the same quadrupole power. 

3. Phase analysis of the WMAP first (WMAP I) and third 
(WMAP III) years the CMB quadrupoles. 

Below we examine the phases of the ILC and those of the foregrounds in more detail. 
Firstly, we plot the phases for the ILC (I) and the foreground and the frequency 



(8) 
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Fig. 1. The map for the ILC(III) quadrupole (top) and the transformed map by the permutation 
operator p (bottom). These two maps have the same quadrupole power. 

maps from WMAPI in Fig. 2. The phases for foregrounds at Q and V band and for 
frequency map at Q-W to show that £2,1 — ^ff (see errors in the Table 1). For K 
and Ka foregrounds the correlation are even stronger than for Q band. 

As one can see from Table 1, the phase of the ILC(I) £2,1 is extremely close 
to the phase of the foreground $2.1 while the phase difference $2,2 — £2.2 is close 
to 7T. To specify the probability of such realization of the phases from uniformly 
distributed and non-correlated CMB phases with the foregrounds, we can use two 
vectors n\ = (cos £2,1, sin £2,1) and n 2 = (cos £2,2; sin £2,2) for the ILC(I) and /1 = 
(cos $2,i, sin <j>2,i) and f% = (cos $2,2, sin $2,2) f° r the foregrounds. 
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As is proposed for the multipole vector in Refs. 22, 26, we can characterize the 
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Fig. 2. Phases of the WMAPI Q, V and W band signals. The phase is represented by the angle 
subtended between the positive x axis and the line. Solid lines arc for (£, m) = (2, 1) modes and 
dotted lines for (£,m) = (2,2). Red lines are the ILC(I) phases and are plotted in all 6 panels 
whereas black lines are for either foreground maps or frequency maps. Notice the overlapping of 
phases at the (2, 1) component between the ILC(I) and Q and V foreground maps, and Q, V and 
W frequency maps. The difference is listed in Table 1. 

correlations of these unit vectors in terms of scalar and vector product: (n*j • /j) and 
(ni x fi)- Following Helling, Shupp and Tesilleanu (Ref. ?), we use the vector product 
for estimation of the probability to get realization with Tlx x /i = sin($2.i — £2,1) as 

^($2,1,6,1) = ic i 2 sin2($ 2)1 (10) 

where — kl/(k — p)l. Since <&2.i — £2,1 *C tt/2, this probability corresponds to 
the phase difference $2,1 — £,2,1 multiplied by a factor of 2. From Table 1 one can 
see that for Q band foreground it is about 5% . 

As is claimed by the WMAP science team, improvement of systematic effects, 
including the gain factor, leads to correction of the foregrounds and the ILC(I) map. 
We show below that this improvement reveals significant correlations between the 
ILC(III), the foregrounds and non-cosmological dipole phases. First, we start from 
analysis of the difference between ILC(III) and ILC(I) map. In Fig. 3 we plot the 
quadrupole map for ILC(III)-ILC(I) in Galactic coordinates. 

Correction of the foregrounds in WMAP III and the LLC (III) in comparison with 
WMAP I leads to the changes of the ILC(III) and the foregrounds phases as shown 
in Fig. 4. In Table 2 we list the phase differences for $2,1 — $2.2 and £2,2 — £2,1- As one 
can see from Fig4, the phases of ILC(III) now are significantly different with respect 
to the ILC(I), while the phases of quadrupole for the foregrounds at (2,2) are still 
shifted by ~ ir. The bottom right panel marked with "DIP" in Fig4 shows the shift 
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Fig. 3. The map for difference between the ILC(III)and ILC(I) (top) and that for difference 
between the K band foreground WMAPIII and WMAPI (bottom). 



of phases for the ILC(III) quadrupole with respect to both non-cosmological and 
cosmological dipole phase ((£,m) = (1, 1)) from ILC(III). 

From K and Ka band of the Table 2 one can write down an empirical equation 

6,i + $2,1 - 6,2 + $2,2 = 4.46963(-0.0147) ; K 
6,i + *2,i 6,2 + $2,2 = 4.42253(+0.0396) ; Ka 

(11) 
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The probability for the phases of the foregrounds and the non-correlated phases 
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Fig. 4. Phases of WMAPIII foregrounds and the ILC(III). The phase is represented by the angle 
subtended between the positive x axis and the line. Solid lines arc for (£, m) = (2, 1) modes and 
dotted lines for (£, m) = (2, 2). Black lines are for foregrounds and red lines are for ILC(III). The 
bottom right panel shows the ILC(III) phases in red, the non-cosmological dipole ((£, m) = (1, 1)) 
phase (black solid line) and the cosmological dipole one (black dotted line). 

to satisfy Eq.(ll) is exactly the modulo of difference between £2,1 + $2,1 and £2,2 + 
$2,2, as shown in the brackets in Eq.(ll). 

Correlation of phases in the form of Eq.(ll) can be easily explained in terms of 
phase correlation after the permutation operation on the phases. Taking account of 
^2.1 = £,2,2 and ^2,2 = £2,1, from Eq.(ll) we obtain ^2,1 — $2,1 — ^2,2 — $2,2 and 
rotation of the c?2,m phases by the angle A = ^2,1 — $2,1 transforms the phases 
^>2,i — A, ^22 — A to the phases of the foregrounds $2,1 > $2.2- 

As is shown in Table 2, the probability for ILC and foreground phases to follow 
Eq.(ll) for uniformly distributed and statistically independent from the foregrounds 
Gaussian CMB signal increase from ~ 0.09 for Q band and reach the maximum 
0.22 for the V band and falls down to 0.12 for the W band. This result seems to 
have natural explanation, since the contribution of the synchrotron emission to the 
Q, V and W bands decrease, and the free-free and dust emission dominate over 
synchrotron emission. 

4. Multipole vectors analysis 

In this section we compare the phase correlation method with multipole vectors 
approach, proposed in Refs. 13, 23. We apply the multipole vectors method for the 
ILC (III) quadrupole and kinetic quadrupole. After this we calculate the multipole 
vectors using the method in Ref. 23. 
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The values of ILC quadrupole ai m are: 

1 m 5R 3 

2 +1.147576980 x 10~ 2 

2 1 -5.329925261 x 10~ 5 +4.864335060 x 10~ 3 
2 2 -1.440715604 x 10~ 2 -1.880360022 x 10~ 2 

After correction of the kinetic quadrupole, we have 

1 m 5R 3 

2 +1.147431601 x 10~ 2 

2 1 -5.358225098 x 10" 5 +4.866961855 x 10" 3 
2 2 -1.440600399 x 10~ 2 -1.880334876 x 10~ 2 

According to Ref. 13 we describe these ai m by multipole vector with coor- 
dinates (-0.562380, 0.815276,0.138036) and (0.970920,0.048491,0.234440) which 
correspond to the points on a sphere with Galactic coordinates (in degrees) 
(I = 124.6, 6 = 7.9) and (I = 2.9, 6 = 13.6) respectively. 

We perform the same conversion for the K band. The input data are 

1 m 5R 3 

2 -2.167431593 x 10+° 

2 1 -1.421663761 x 10" 1 4.126564786 x 10~ 2 
2 2 +4.031517506 x 10" 1 1.558263451 x 10" 1 

After correction for the kinetic quadrupole, we have 

1 m 5R 3 

2 -2.167433023 x 10+° 

2 1 -1.421666592 x 10" 1 4.126827419 x 10~ 2 
2 2 +4.031529129 x 10" 1 1.558265984 x 10" 1 

This gives us the corresponding multipole vector coordinates and the Galactic 
positions on a sphere: (-0.560275,0.076989,0.824721), (I = 172.2,6 = 55.6) and 
(-0.479820,0.113498,-0.869995), (/ = 166.7,6= -60.5). 

Using statistics for such vectors, we have found area vectors in a form 

U> (2;1,2) = „(2,1) xu (2,2)_ {u) 

This gives us the area vectors with coordinates Win: (0.1850,0.2657,-0.2194) or 
{I = 55.2, 6 = -34.1) and w K : (-0.1602, -0.8826, 0.0033) or (I = 259.7, b = 0.2). 

After this we calculated statistics following Ref. 23 to estimate correlation prop- 
erties: dot product of ILC and K-channel vectors 

= vf/p • v ™ , (13) 

dot product of ILC and K-channel area vectors 

(2:1,2) (2:1,2) 

A(2;l,2) _ W III ■ W K (u) 
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and ratio of lengths 



|7/, (2;1 ' 2) 
(2;1,2) _ \ W K 



(15) 




So, we have got 



d^ 2 ' 1 ) +0.4908 

d< 2 ' 2 ) -0.6643 
A (2;i,2) „ . 7550 

r (2;l,2) +2 .2939 



As one can see from this Table, all the dot and cross-products are not especially 
close to unity (see for comparison Table 2 in Rcf. 23 for ILC(I) and ILC(III), for 



which S- 2 -^ ~ 0.973,d( 2 ' 2 ) ~ 0.956, A< 2;1 < 2 ) ~ 0.955 and r( 2;1 < 2 ) ~ 0.952), which 



means that the multipole vector approach is not very sensitive for analysis of the 
CMB-foreground correlations. 

5. Discussions 

5.1. Systematic effects ? 

For interpretation of correlations between the WMAP low multipoles (i.e. alignment 
for £ — 2,3, planarity for i = 5 etc.) we need to know how possibly systematic effects 
determine their properties. This problem will be even more crucial if one is to explain 
the deficit of the ILC quadrupole power. We will need to implement modification 
of the theory of inflation (cut-off in power for primordial adiabatic fluctuations at 
spatial scales > 10 3 Mpc), the Bianchi Vllh model, or any theoretical changes of 
the cosmological scenario, assuming primordial origin of the ILC(III) low amplitude 
quadrupole. 

Probably, the best and more significant illustration of the systematic effects is the 
improvement of the gain factor performed by the WMAP science team, which could 
be the reason why the quadrupole from ILC (I) has different amplitude and phases 
in comparison with from ILC (III). More importantly, the quadrupole power of the 
ILC (I) and the ILC (III) does not change as much as the phases reveal of a new sort 
of correlation between the ILC (III) phases and the foreground phases (see Eq.(ll)). 
In Fig. 5 we plot the corresponding phases for the map of difference between the 
ILC(III) and the ILC(I) quadrupoles. Next type of systematic effect is illustrated 
by the comparison of the ILC (III) quadrupole and the quadrupole derived by de 
Oliveira-Costa and Tegmark (Ref. 35) from the WMAP I and the WMAP III. Differ- 
ence between the ILC(III) and de Oliveira-Costa & Tegmark (Ref. 35) quadrupoles 
is mainly related with different Galactic mask, which we will call M05 and M06 
(following Ref. 35). As one can see from Fig. 5 and Table 3, the error of the £2,1 
ILC phase reconstruction is about 0.1 radians for M06 and it is negligible for M05 
mask. For £2,1 the corresponding error is <~ 0.12 radians (M06) and ~ 0.05 radians 
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Fig. 5. Phases of WMAPIII foregrounds and the ILC(III). The phase is represented by the angle 
subtended between the positive x axis and the line. Solid lines are for (£, m) = (2, 1) modes and 
dotted lines for (I, m) = (2, 2). Red lines are for ILC(III), black and blue lines for K and W band 
foregrounds, respectively. Green for ILC(III)-ILC(I). Thick red dash line is for non-cosmological 
dipole. Thick blue and green are for M05 and M06 respectively, in which dash lines and triple 
dot-dash lines are for (2, 1) and (2, 2), respectively. 





ILC(III) 


A 


M05 


M06 


(2,1) 


1.582 


5.856 


1.528 


1.668 


(2,2) 


4.059 


4.909 


4.010 


4.181 



(M05). This means that with uncertainties about ±0.12 radians Eq.(ll) is correct 
for all the K-W foregrounds. 

Let us discuss another important correlations of the ILC(III) phases and the 
phases of non-cosmological and "cosmological" dipoles shown in Fig. 4. By dipole 
phase we mean the phase of (£, m) = (1, 1) mode. The angle between £2,1 = 1.528 
and the non-cosmological dipole \& nc — 4.824 is about a — f nc — £2,1 = 3.295 
radians, while the angle between £2,2 = 4.0586 and the phase of cosmological dipole 
* c = 0.8379 is = £2,2 - * c = 3.2207. Thus, a - ~ 0.074 rad for the ILC(III). 

We would like to point out that the above-mentioned correlations between the 
phases of the ILC(III) and those of the WMAPIII foregrounds are mainly related 
to the changes of the gain factor. The crucial part of these correlations is the angle 
between £2.2 — £2,1 an d for the foreground phases <E>2,2 — $2.1- Note that transition 
from the WMAP I to WMAP III leads to rotation of the phases for the ILC and 
for the foregrounds as well. Since WMAP III data seem more accurate in terms of 
systematic effect removal, decreasing of direct correlations between £2,1 and $2.1, 
typical for the WMAP I, reveals another type of correlation between the ILC(III) 
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and the foregrounds. In particular, in addition to Eq.(ll) one can get 

^2-^1-1-186 rad. 

&,i ~ £a - 1-133 rad. (16) 

which means that the phases for the ILC(I) and ILC(III) have some regular changes. 
It is quite possible that these changes reflect some correlations between the LLC (I) 
phases and the WMAP I foregrounds. The comparison between synchrotron phases 
for the K band shows that from Eq.(16) one can obtain ^ 2 — £2*1 — 2.31907 radians, 
while $21 — $2,2 — 2.34936 radians. So, with accuracy about 0.030 radian we 
get : £"1 — C2.2 ~ ($2,1 — ^2,2)- C* n the other hand, from Eq.(16) one can get 
£™ — 2^2,1 — £2,2- Finally combining these two equations we get: 

2^,1 + $2,1 - 2^2,2 + $2,2 (17) 

As for the WMAP III, one can see that the phases of the ILC(I) and the WMAP I 
foregrounds per K band are strongly coupled in the form of Eq.(17). However, note 
that the Eq.(17), like the Eq.(ll), is "empirical", displaying possible relationship 
between the phases of the ILC and the foregrounds, which requires explanation in 
the framework of the CMB and foreground separation scheme. 

Implemented in this section phase approach can be reformulate in terms of stan- 
dard correlation between different components of quadrupoles. Since the quadrupolc 
component of the signals (the CMB and the foregrounds) contain only 02,1 and 02,2, 
the coefficient of the cross-correlation between them is simply 

w ilc |a2,i||a2,2|cos(&,i -62,2) , . , . 

K l,2 = /, , 9 | , f = COS(^2,l - 42,2) (18) 

vl a 2,ir|a2,2r 

Even if the ILC quadrupolc is of Gaussian nature, for a single realization of the 
random process, as our Universe does, this cross-correlation coefficients are not 
informative apart from some peculiar values such as or 1. However, combining the 
coefficient with the same coefficient K( 2 , defined for the WMAP foregrounds 

by simple comparison between these two numbers we can conclude about common 
morphology of these signals. 



5.2. Primordial origin ? 

As is mentioned in Introduction, starting from the COBE experiment, the proper- 
ties of the CMB quadrupolc and particularly, the deficit of its power in comparison 
with the WMAP best-fit ACDM cosmological model has attraction of serious at- 
tention. Described in the previous section, phase correlation between the ILC and 
foregrounds contain important information beyond the power spectrum, providing 
significant restrictions on modification of the theories of inflation. From our anal- 
ysis it is clear that models of inflation based on primordial Gaussian fluctuations 
of the inflaton field can not explain the phase correlations. Moreover, the question 
is if these cross-correlations can be explained in the framework of widely discussed 
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the Bianchi Vllh anisotropic cosmological model (see, e.g. Rels. 24, 39 25) as this 
model predicts specific properties of the phases for the CMB quadrupole ? To an- 
swer this question we extract the quadrupoles of the maps derived in Refs. 39 and 
25. In Fig. 6 we plot both CMB maps for the range of multipoles I < 10 and their 
corresponding quadrupole components. It can be seen clearly that morphology of 
these quadrupoles is different from that of the ILC(III) quadrupole (see Fig.l). 

To characterize this difference, in Fig. 7 we plot the phase diagram for the 
quadrupole modes of the ILC(III) and those from the above maps. 

In spite of differences in phases between the ILC(III) and both Bianchi Vllh 
quadrupoles, they have remarkable correlations. Namely, the phase difference for 
the ILC(III) phases is Am = £2,2 — £2.1 — 2.477 rad, while for Ref. 39 quadrupole 
it is A j = £2 2 — £jf.i — 2.537 radians and for Bridges et al. (2006) quadrupole it 
is Ab = £22 ~ £2^1 — 2.560 radians. As one can see, with error about 0.07 radians 
these angles Am, Aj and Ab are the same, which means that by simple rotation of 
the phases we can get the same morphology of the maps. In particular, by rotation 
of the Ref. 39 quadrupole phases by the angle f3j = 0.367 radians counter-clockwise 
and by rotation of Bridges et al. (Ref. 25) phases by the angle (3b — 0.317 radians 
clockwise we get the LLC (III) phases with above-mentioned accuracy. 

In Fig. 8 we plot the difference between the Bridges et al. quadrupole (Ref. 25) 
after rotation of phase by the angle /3b = 0.317 radians and the ILC(III) quadrupole. 
Note that the peak to peak amplitude of the maps shown in Fig. 8 are at the same 
range as the ILC(III) quadrupole: at —200, 200/iK . That means that in spite of 
correction of the phases for Bridges et al. quadrupole (Ref. 25), the amplitudes of 
fl2,o -i- a-2,2 modes still are not optimal in comparison with the ILC(III) amplitudes. 

6. Conclusion 

We have presented the phase analysis of the LLC (I), ILC(III) and de Oliveira-Costa 
and Tegmark (Ref. 35) quadrupoles in comparison with the foreground phases. 
We have shown that the ILC(III) quadrupole has strong correlation with the fore- 
ground phases, mainly with the synchrotron emission per K and KA band fore- 
grounds. We have checked out the possibility that the WMAP third year data 
release quadrupole can be explained by implementation of the Biachi Vllh model 
discussed in Refs. 39 and 25. All these models need additional corrections of the 
phases for a 2 ,i and a 2 ,2 components in order to math the phases of the ILC(III). By 
analysis of the de Oliveira-Costa and Tegmark (Ref. 35) quadrupoles related with a 
different Galactic mask we can conclude that detected above correlations between 
the ILC and foregrounds weakly depend on the type of the mask, if the area covered 
by them is small in comparison with the 4tt. We have shown that transition from 
the ILC (I) quadrupole to the ILC (III) changes significantly the properties of the 
CMB-foreground phase correlations through renormalization of the gain factor. We 
believe that this is the main reason for discussed above phase correlation, clearly 
demonstrated importance of accurate removal of the systematic effects. 
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Fig. 7. Phases of the quadrupole modes of the ILC(III) and the Bianchi Vllh model. The phase 
is represented by the angle subtended between the positive x axis and the line. Solid lines are for 
(£, m) = (2,1) modes and dotted lines for (2,2). Red lines are for ILC(III), black for Jaffe et al. 
and green for Bridges et al. 
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